AD-AlQl 999 MASSACHUSETTS INST OF TECH CAMBRIDGE DEPT OF OCEAN E—ETC F/G 20/11 
CRUSHING ANALYSIS OF ROTATIONALLY SYMMETRIC PLASTIC SHELLS.<U> 

JUN 61 J G DE OLIVEIRA* T WIERZBICKI N00014-80-C-0616 

UNCLASSIFIED 81-8 NL 




























6 6 6 101V QV yjy 4ng 





* On leave from the Institute of Fundamental 



6 

81 7 24 020 







-i- 

INDEX 

Abstract . 

1. Introduction . 

2. Analysis of Discontinuities . 

3. Mechanisms of Plastic Deformation . 

4. Determination of the Crushing Force .... 

5. Tube Inversion . 

6. Crushing of Conical Shell . 

7. Spherical Shell Under Point Load . 

8. Crushing of a Spherical Shell Between Rigid 

Plates . 

9. Rigid Boss Loading . 

10. Spherical Cap Under External Pressure Loading 

11. Effect of End Condition . 

12. Conclusions . 

Acknowledgement . 

References . 

List of Figures . 


Figures . 

Library Card . 

Report Documentation Page 


Accession For 

•NTiS • CR'SI * • !*/ 

OTIC T'S r- 

4J;iennou>. .■ ">d • • 

Just if 1cetlcn_ 


By-__- 

Distribution/ 
Availability Codes 
Avail and/or 
'Diet Tpectal 


ft 






























- 11 - 


Abstract 

_ J 

The crushing analysis of rotationally symmetric plastic 
shells undergoing very large deflections is presented. A 
general methodology is developed and simple closed form solu¬ 
tions which can be useful for practical applications are derived 
for the case of a conical shell and a spherical shell under 
point load, a spherical shell crushed between rigid plates and 
under boss loading, and a spherical cap under external uniform 
pressure. The effect of the end conditions and the limitations 
of this approach are discussed in detail. 
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1. Introduction 

Thin metal shells of revolution subjected to axisyiranetric 
compressive forces are capable of carrying substantial loads 
when deflections exceed by one or two orders of magnitude the 
thickness of the shell, and become comparable to the largest 
linear dimension of the structure. Thus, the deformation 
process is far beyond what is commonly understood as post- 
buckling behavior of elastic or plastic shells. A distinctive 
feature of such deformation mechanism is that the strain energy 
function (in the case of elastic shells), or dissipation function 
(in the case of plastic shells), is concentrated in narrow zones 
(fold lines or hinge lines) while the remainder of the struc¬ 
ture is undergoing a rigid body motion. Furthermore, with some de¬ 
cree of symmetry, the fold lines are forced to move as the deforma¬ 
tion process goes on. Depending on the geometrical parameters 
of the structure and imperfections, either a symmetric or un- 
symmetric deformation pattern may develop, even though the loading 
and the structure itself are axially symmetric. 

A known example of non-symmetric but regular deflection field is 
the Yoshinura buckle pattern in conpressed thin elastic shells. In 
the range of very large deflections, the geometry of such a field 
can be conveniently described using the concept of isometric 
transformation of surfaces. The foundations of the respective 
mathematical theory were laid by Pogorielov, [1], who also 
presented brilliant applications of his theory to the analysis 







of post-critical behavior of thin elastic shells, [2], His 
work was later followed by Lukasiewicz and Szyszkowski [3], 
while Foster [4] presented a parallel but independent study 
of the compression of thin tubes. Pugsley and Macaulay 
appear to be the first to study the crumpling process of 
thin cylindrical shells with diamond shape lobes in the plastic 
range [5]. Although no provision was made to accomodate in 
the theory the travelling hinges, the corresponding calcula¬ 
tions for the mean crushing force, especially the more recent 
ones [6], agree well with experimental results. Further contri¬ 
bution to the understanding of the process of progressive 
crumpling of tubes in the multi-lobe modes was made in [7] 
and [8], but no satisfactory solution of this interesting problem 
was offered to date. 

As tubes get thicker (R/h < 50) , a well-known transition 
takes place from the diamond pattern to the crinkling (concertina) 
mode of deformation, discussed by Alexander [7], Assuming a 
simple deformation mechanism with circumferential extension 
and equating the rate of internal energy dissipation to the 
rate of work of external forces, he was able to calculate an 
average crushing force. Minimizing next the force level he 
evaluated the length of the local buckling wave, and derived a 
simple expression for the optimal value of the crushing force. 

The above procedure has proved to be most effective in the 
approximate analysis of the crumpling process of various thin- 
walled structures [8,10]. 










The tube inversion (outside-in or inside-out) is another 
example of an axially-symmetric deformation of tubes. Two 
solutions for the magnitude of axial loading necessary to main¬ 
tain the plastic flow are available. One, due to Al Hassani 
et al. [10] , is good as long as the assumption of a constant 
thickness holds, and another one, due to Abraraowicz [ 11 ], takes 
into account changes of thickness and circumferential curvature. 
The problem of hemispherical shell or spherical cap loaded by 
central point-load, rigid boss, or crushed between plates, received 
much less attention in the literature than the related buckling 
problem under external pressure loading. Updike and Kalnins 
performed a thorough analysis of postbuckling behaviour of 
elastic spherical shells compressed between rigid plates, [12, 

13] . Using F.eissner equations for shallow rotationally symmetric 
shells, they determined the force-deflection relationship in 
subsequent stages of deformation, and the point of bifurcation 
into a symmetric inward dimple, and later into an unsymmetric 
shape. In the subsequent study [14] Updike presented an approxi¬ 
mate solution for the rigid-plastic shallow shell undergoing 
symmetric deformations. An extensive theoretical and experi¬ 
mental investigation into this problem was made by Kitching 
et al. in [15], where the transition into an unsymmetric plastic 
mode was calculated numerically. Experiments in which polygon 
rather than circular circumference dimple were observed, were 
reportedly made earlier by Pian [9]. 









Leckie and Penny [16] performed a series of tests on 
carefully manufactured hemispherical shells loaded centrally 
by a rigid boss. These experiments were followed by a theoretical 
study due to Morris and Calladine [17] , in which the deformed 
shape of the shell and the force-displacement characteristics 
were obtained by a sequence of upper bound calculations , in 
conjunction with an ingeneous idea of treating the shell as a 
three-dimensional body, [18]. 

The work of Morris and Calladine consititutes a milestone 
in the understanding of the crushing behaviour of shells of 
revolution. It shows through relatively simple numerical cal¬ 
culations that plastic deformations are indeed confined to a 
relatively narrow ring or section of a toroidal surface, and 
that the ring is moving outward as the deformation process 
goes on, leaving a rigid region behind. A snap-through effect 
was observed whose extent depends on the boss size p. 

The plastic behaviour of shallow spherical shells under 
uniform pressure loading is similar to that of a shell under 
rigid boss loading, as far as the occurrence of the plastic 
snap-through phenomenon is concerned. This problem was studied 
by Puszek [19] and Jones and Ich [20], using an approximate yield 
condition for rotationally symmetric shells. Recently, the 
same problem was reconsidered by Kondo and Pian [21] by applyine 
the generalized yield line method, v/hich takes into account 
the changes in the geometry of the structure. 









In the present paper we shall extend the approach developed 
in [7] and [17] to general shells of revolution, with the 
purpose of deriving some simple results useful for practical 
purposes. A detailed discussion will be presented on the 
restrictions imposed on the solution by the conditions of 
kinematic continuity on the moving hinge circles, and inextensi 
bility of the shell in the meridional direction.. Introducing 
simple displacement and velocity fields, closed-form solutions 
on the crushing force will be derived for a variety of struc¬ 
tures, including conical and hemispherical shells compressed 
between rigid plates, and a spherical shell under a point-load 
or central boss loading. The influence of the boss size, 

and clamped edge conditions will also be discussed at 
length. 

The analysis will be restricted to axisymmetric deforma¬ 
tion modes only. A simplified approach to study the lateral 
crushing of cylindrical shells using a polygon shaped fan of 
flat elements was developed by Morris and Calladine [22]. 

This method is particularly suited for treating nonsymmetric 
but regular deformation modes in hemispherical shells such 
as the ones observed by Pian [ 8 ]. The treatment of non¬ 
symmetric problems with smooth deformation fields will be 
the subject of a future publication [23]. 








2. Analysis of Discontinuities 

The condition of kinematic continuity on the moving hinges 
or hinge circles in rigid-plastic beams and flat circular plates 
were formulated by Hopkins [24] and Hopkins and Prager [25]. 

Before writing the continuity equation for a general shell 
consider first the one-dimensional case, since the physical 
arguments leading to these equations are the same. 

Let T denote a straight line yield hinge whose position 
on the flat surface element is described by a scalar equation 
x = v(t). The downward deflection and velocity of the material 
point lying on the shell middle surface are denoted respectively 
by f(x,t) and f(x,t). We assume that deflections and velocities 
are continous at any time and at any point of the shell. 

Since the deflection f is continuous across F, the time deriva¬ 
tive of f(v(t),t) along T is also continuous across F. Thus 

[f] + 0 (t) [f' ] = o (i) 

where the brackets denote discontinuties, if any, of the enclosed 
quantities across r, i.e. [n] - n + - n_. Now, f is continuous across 
T and Eq. (1) therefore shows that the slope f' can be discontinuous 
only across stationary yield hinges, for which v(t) =0. If the 
yield hinge T is not stationary, then f’ must be continuous across 
T. By repeating similar arguments one can show that 

[f'l + v(t) [f”] = 0 (2) 

provided v (t) ^ 0. 







Consider now a hinge line moving with a constant velocity 


v (t) = V down the undeformed flat surface element 

leaving a plastically deformed region, Fig. 1. The rate of 

rotation and curvatures on both sides of the moving hinge are 

respectively f + = 0, £'_ = 0, f+ = k, f" = 0. Now, Eq. (2) 
reduces to 

0 + VK = 0 (3) 

tve can see that the moving hinge line imposes a constant 
curvature < = 4 on the initially flat surface element if the 
slope of the velocity field 6 is related to the hinge velocity 
by (3T\ 

Return now to the general case, following the derivation 
presented in [11] and consider a curved discontinuity line F 
moving down the middle surface of the shell x* The line r 
is dividing the surface x into two separable parts 
X = ( + ) XUrU ( ^x* Then ^X denotes the deformed part, subjected 
to plastic deformation. The vector field f is defined over the 
entire surface x- Now, the previous analysis still applies 
except that the first and second derivation of f in (1) and 
(2) should be replaced by the evalutation of the first and 
second gradient of f on the surface x. This is equivalent 
to the replacement of partial derivatives in (1) and (2) hy 
covariant derivatives, so these eauations now read 


Tf 1 ! + v f 

fi | 1 

L J n [ 

. a J 


(4) 


0 


i = U,2,3} 
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where n is a surface vector normal to the discontinuity line 
r, while V n denotes the velocity of the hinge line in the n 
direction. 

Considerable simplification can he obtained in the case 
of rotationally symmetric shells. The hinge line becomes now 
a circle and each of the continuity condition (4) and (5) 
yields only one scalar equation in the meridional plane 

[f] + v ff, ] = 0 (C) 

n a 


[f, a ] 


v [f. 

n act 


where [f, ] is a jump in the rate of rotation and [f, 1 is 

a -it- aa 

the corresponding jump in the principle (meridional) curvature. 

Take now an intermediate stage of the crushing process 
’•-'here plastic deformations are confined to a relatively narrow 
zone contained between an inner and outer hinge circle B and C, Fig. 
The outer hinge circle B moves down the undeformed shell with a 
velocity V imposing a larger curvature, i.e. changing the 
local radius of curvature from R_ to r F . The material in front 

d b 

of the hinge is rigid and undergoes a constant and finite rate or. 
rotation 0. Thus, the continuity condition reduces in this 


case to 

~“ b + v 44 ' h) = 


The inner hinge C is travelling with the velocity and remover 

the curvature, or strickly speaking, it changes the small radius 

of curvature r back to the large one K . be can then write 
C C 








a similar equation to (8) expressing the continuity constraint 

u + V I — - = 0 ( 

c c( r c b c ; 

For the material behind the circle C to be rigid, the shell 
should translate rather than rotate. Hence, there must be at 
C a counterrotation ^ of the same magnitude as at B. The 
condition ojg=o c "‘implies the following relation between the six 
kinematic and geometric parameters of the problem 



Equation (10) can be further simplified. If for example the 

deflected shape of the shell is assumed to be isometric to the 

original one, it must be its mirror reflection i.e. Rg = R c 

see [1]. Furthermore the shape of the shell middle 

surface between the hinge circle B and C may be approximated by 

a section of a circular toroidal surface so that r_ = r = r, 

B c 

as suggested in [14,11]. 

Equation (6) carries the information that the shape of shell 
at B and C should be continuous at all times, except at the 
stationary hinge, i.e. at the rigid boss or clamped edge. 

In the next section we shall see whether conditions (6) 
and (10) can be made compatible with other simplifying 
assumptions regarding the kinematics of the crushing process. 

3. Mechanisms of Plastic Deformation 


The foregoing analysis has shown that there are two con¬ 
centrated plastic hinges at B and C and that the cross section 
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be tween them rotates as a rigid shape with an angular velocity w, 
as the displacement of the shell increases. Our aim is to evaluate 
the internal energy dissipated in such a. mode of deformation. 
Equating this energy to the rate of work of external forces 
an instantaneous crushing force can be found necessary to 
maintain the plastic flow. 

The dissipation is due to a discontinuous velocity field 
at plastic hinges, and to a continuous deformation field in the 
shell section between B and C. Consider first the latter.. 

In rotationally symmetric shells there are four components 
of the generalized strain rate and stress fields. The meridional 

curvature rate < is infinite at hinge circles and vanishes 

•D 

between them. This follows from the assumption that the rate 
of rotation is a step function. 

It should be noted that the meridional curvature itself 
< 5 , governed by (8) is kept constant at any material point 
passing through the plastically deforming annulus, but does net 


have to be constant along spatial coordinate between B and C. 


The circumferential curvature rate < e is continuously 

changing between B a nd C from positive — through zero to negative 
i_ r ' B 

R r/ and consequently some energy is dissipated in this deformation 


mechanism. An exact formula for k p was derived in [11] using 
a rigorous Eulerian description. In the same paper it was 


shown that the associated energy dissipation is a small fraction 


of the total energy. 
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For the sake of simplicity we shall disregard < in the 

9 

equation for the energy dissipation- A very good estimate of 
the contribution of ic to that energy can be given in a 
straightforward manner. 

Suppose the circumferential curvature is reversed between 
B and C from 4 to . The rate of work done is twice 
the energy required to remove the curvature from 4 to zero. 

The curvature is removed by rotating two ends of a cut circle 
through the angle 2rr, (cf. [26]). Thus, the rate of energy 
in the considered mechanism of deformation can be approximated 
by : 

E(k 0 ) = 2 (2tt) M q V = 4 ttM q | w (11) 

where it was assumed that there are no interaction between the cir¬ 
cumferential bending moment and the remaining non-vanishing 
components of the generalized stress field. The velocity V in 

(11) is an average downward velocity of the deforming zone, 
b • 

V = j to; see Fig. 3. 

At certain stages of the crushing process considerable 
meridional membrane forces may be developed. In the present 
paper the meridional extension rate is assumed to vanish in 
the zone of continuous deformations. Any extension or 
shortening of the arc length in the meridional direction is 
accommodated by the hinge circle B and C. Thus, we do not 
require inextensibility of the material in the meridional 
direction but instead extensional plastic hinges are intro¬ 
duced. This assumption is in the spirit of Calladine's: approxi¬ 
mation in the problem of moderately large deflection of clamped 
circular plates [18]. 
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The only significant component of the generalized strain 
rate is the circumferential compression A . In order to calcu¬ 
late the associate dissipation energy v;e shall strictly follow 
the approach developed in [17], and treat the shell section 
between B and C as a three-dimensional tody. The rate of internal 
energy dissipation is equal to 


E = 


/Ve dv 

V 



a 

o 


x e 9 dS 


where V and S denote respectively the total volume and area 
of the cross section between B and C, and x denotes- a current 
distance of the given point F from the axis of symmetry of the 
shell, see Fig. 3. Note, that the yield stress is positive in 
tension zones and negative in compression zones. The circumfer- 
ential strain rate is given by Eg = where y is the current 
distance from the axis of instantaneous rotation I-I. The 
position of this axis with respect to the mell cross section 
is yet to be determined. Substituting e a in (12) yields 


E. . 
int 


= 2710 


'o-/l 


y |cS 


(13) 


One has to distinguish now two special caeec. If the rise of 
the toroidal surface over the points L anc. C is comparable 
to the shell thickness, the assumption about non-uniform 
distribution of t a across the thickness is justified, and (13) 
holds. If the converse is true, one can replace dS = hdx, and 
(13) simplifies: 
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F = 2-rro hu> /|y!dx (14) 

iht ° bc 


The integral appearing in (14) is the area between the 
middle surface of the annular section of the shell and the 
axis of instantaneous rotation, as noted in [14]. 

The energy dissipated in any of the two plastic hinges 


is 


h 


h 


• disc 
E int 


= 27ia 


/ a o i 0 dy = 2 ira U 0 o / |y|dy 


( 11 ) 


2 2 

where a denotes the position of the hinge with respect to the 

axis of symmetry, and u@ is in-plane component of the velocity field. 

One can easily evaluate (15) as a function of the position of 

the center of instantaneous rotation with respect to the middle 

surface of the hinge £: 


£ ^ for £ = 0 (16a) 

£ 2 + | 2 > u < £ < | (16b) 

£h 1 for £ > ^ (16c) 

In the present approach the hinge circle is indeed a 

generalized plastic hinge analysed for example 

by Jones [27] , since either bending or extension or both can 

be devleoped in it. 



The total dissipated energy is the sum of energies 
dissipated in continuous and discontinuous dissipation fields. 
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4. Determination of the Crushing Force 

In the actual solution the rate of change of external work 
is always equal to the internal dissipation energy. A similar 
equality for any kinematically admissible velocity and dis¬ 
placement field can be understood as a definition of the approx¬ 
imate value of the external loading. The more close the 
assumed fields are to the real one, the better the approximate 
solution is. A considerable success of this m.ethcd in predicting 
the crush behavior of various thin-walled structures (see [7,15]) 
is due to the correct choice of the admissible field with few 
free parameters, and also to a somehow weak dependence of the 
solution on the details of the assumed share of the shell. 

For example, any shape of the shell middle surface between 
inner and outer hinge circles, satisfying kinematic continuity 
and giving the same area in (14), would yield the same force 
level. In the present paper we shall take the shape of the shell 
middle section in the form of a circular or parabolic arc. 
Furthermore, we shall assume the deformed but rigid portion CD 
to be a shifted mirror reflection of the original shape AB 
(doted line on Fig. 2) by a magnitude b. This implies opposite 
but equal in magnitude radius of curvature and slopes at points 
B and C. From the continuity of displacements and slopes at 
B and C, it follows that the transition shape is determined to within 
a single parameter x>/hich in either the distance between points 
B and C, or local radius of curvature r. Another free parameter 
is the position of the axis of rotation 
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The rate of work of external loading is 

13 * 

P -k— uj - for crushing between plates (17a) 

Pbw - for point load or boss loading (17b) 

Now for any given value of the parameter describing the 
crushing distance (for example central deflection or central 
radius of the annular zone), the integral (14) is evaluated and 

the value of the force P becomes a function of two free para¬ 

meters £ and b. 

The function P( b,£) was studied and it was found that for a 
certain value of t, the force attains a minimum. Cn the other 
hand, except at a very early stage of the crushing process, P 
attains a lower limit rather than a minimum with respect to 
£ = 0, Fig. 5. 

Assuming the minimum principle to hold for geometrically 
nonlinear problems, the lowest possible value of P is identified 
as the best approximate solution. This approach is equivalent 
to the sequence of upper bound numerical calculations performed 
by Morris and Calladine [17] . 

5. Tube Inversion 

The approximate solution of the stationary plastic process 
of tube inversion with various degrees of approximation can be 
found in [11] and [15] . We shall study the same problem to explain 
the details of the derivation and to determine the accuracy of 
the present approach by comparing it with known results. 
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Consider a tube of the initial radius R being inverted 
inside out to a final radius R + 2r. At the hinge circle B 
the straight line generator of the cylinder is bent into the 
curvature p and at the hinge circle C the curvature is re¬ 
moved, Fig. 6. A simple and yet accurate solution of this 
problem, accounting for continuous circumferential extension and 
concentrated bending at B and C, was worked out by Kitching et al. [15] 
and has the form 


P min . 1 3r ~ , QO fR. fP-h 

2rtT * 4 Vt ’ 6-92 Vh' r "VT < 1S) 

In order to apply the present approach, we assume that the 
portion of the shell between hinge circles is and C form one 
naif of the toroidal surface. We also assure for simplicity 
chat r/h > 1, so that Eq. (14) applies and t, -• 0. The distance 
between B and c is h = 2r and the integral appearing in (16) 
represents the area of the half circle. Using this result and 
Eq. (15), the rate of energy balance equation yields 

•2 ’ . h 2 

Pr oi = 2 t to hw—s— + 2 tto_uA- [R + (R+2r) ] (K ) 

O 2. 0 4 

The rate of infinitesimal rotation w drops out from both sides 
of (19) and the dimensionless force becomes 



(?-C) 


f’inimizing the right hand side of (20) with respect tc r and 
substituting the optimal radius back into (2C) we can obtain 


the final formula: 
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j2i2 - 4^2 + 2 = 7.09^ ♦ 2.- r =Vt ,21 > 

The above results compares well with solution (18), especially 
for large R/h. 

It is interesting to estimate the effect of the variable 
circumferential curvature, according to the approximate formula 
(11). One can easily find that the magnitude of the non- 
dimensional force in (21) increases then by 2, which again is a 
small contribution for large R/h. At the same time the optimum 
radius of inversion remains the same. 

It can be checked that the condition (10) of kinematic 

V V 
B C 

continuity, reduces now to — = — and thus is satisfied. 

Finally it should be noted that the crushing force of a 
tube loaded by a central force rather than crushed between 
plates, is one half of that predicted by (21). 

6. Crushing of a conical shell 

The geometry of the partially crushed conical shell is shc'x 

in Fig. 7. We assume that the portion of the shell subjected to 

localized plastic deformation forms a section of a toroidal c urface 

and that the inverted cone is of the same angle a. The area of 

r 2 

the section of a circle between the points E ; and C is S = y [(jr-2a) 
- sin(7r-2a)] and the distance b = 2rcosa. It is now straight¬ 
forward to work out the solution which takes the form 

P — 

Vh' V (7T " 2a) ' sin(7T - 20l) (22) 
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-V 


ha 


(iT-2a) - sin(ir-2a) (23) 

Here a is the measure of the crushing distance, which is related 
to the central deflection by 


w = 


2a 

tana 


(24) 


The formulas become increasingly inaccurate for a -r- and 
also for small since then the assumptions leading to 5 = 0 
are not satisfied. In these cases a more exact solution can 
be obtained when needed, using (13) rather than (14) and 
minimizing the solution with respect to 

Equation (23) predicts an increasing radius of the toroidal 
surface, r > 0, as the crushing zone is expanding. The change of 
the radius is related to velocities of hinge circles by 


2r = (V B ~V c )tan a 


(2i) 


At the same time we have assumed that = R c = » and r g = r^ 
which yields V = V . We see that the condition of kinematic 

D L. 

continuity is not satisfied by the present solution. The 
difficulty encountered can be overcome by considering the general 
case in which r_ f r_ R_ ^ ®. Now, the toroidal surface ceases 

D v- C 

to be of a circular cross-section and the hinge circles B and C 
will no longer be positioned on the same level. As a result 
membrane force will develop in the hinge C which would result in 
a slight increase of the force level. Also, with a variable 
radius R (a) the hinge circle will produce a curved line rather than a 
conical one on the meridional olane. It is not difficult to write 
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a system of ordinary differential equations with a shifted 
argument describing variation of V n V_ r_ r_ and R_ with 
the crush distance. We shall not attempt to study this sytem 
since it is believed that insiginificant changes would have then 
been introduced to the solution. The problem however is 
interesting by its own right. 
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Equations (26) and (29) furnish a parametric representation of 

the function force versus displacement. A plot of this function 

for several values of the radius to thickness ration R/h is shown 

in Fig. 9. The dependence of the solution on R/h appears to be 

h 

weak. Indeed, treating ^ as a small parameter and expanding the square 

h h 

root in (29) in power series with respect to — arcund — = 0 one gets 

2 2 

1 - 4 ! —- 1 - 2 R + • • • (30) 

so that (29) can now be approximated by 

_ . 2 

w ^ sin ct /n1N 

h 2 K a-cosa) - 2 —(31) 

Combining next (26) and (31) to eliminate the parameter R/h 
we obtain 

w _ / p min \ 1 1 -ccso -l n ( sin ? 'o ) 
h \2 M o tt + 7 |2n(a)j z | n (a) j 

Expanding the trigonometric function appearing in (32) in power 1 
series and retaining terms up to the third power, the two 
coefficients in (32) become 


1-cosg _ _1 2sina - sin2g 
2n (a) 4 2a - sin2a 


2sin 3 a 2s 


in 3 a a 3 (, a A 
- sin2a 2 \ 6 / 


(33a) 


(33b) 


Thus, the first coefficient turns to be constant, while the 
second varies very little which a, especially for a < 1. Sub¬ 


stituting (33) back into (32) , the final formula for the force- 
deflection relationship takes the form 










The plot of the function (34) is given in Fig. 9 by a broken 
line. The agreement with the exact solution (full line), which 
depends on R/h, is seen to be very good. 

It is interesting enough that the same solution can be 
obtained by making a different set of assumptions. Suppose the 
hemispherical shell is approximated by a parabolical shell. The 
relationship between the angle a and central deflection is now 
given by 


w 

h 


R 

h 


2 


sin a 



(35) 


Replacing likewise the circular cross-section of the toroidal 
surface by a parabolic arc 



the area integral in (14) can be readily evaluated, and the 
expression for the crushing force takes the form: 


(36) 


2^rM, 


2b 2R sina i 

- 3 h Slna + -E- 1 


(37) 


Optimization the right hand side of (37) with respect to b 
’-esults in 


b = 3hR 


(38) 


Now combining (35), (37) and (38), the expression (34) for the 

crushing force is obtained in the same form as before. This result is 
important from the practical point of view, since the mathematics 
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involved in the alternative derivation of (34) are much simpler. 
We shall take advantage of this property, and in the remainder of 
the paper dealing with spherical shells always replace a circle 
by a parabolic arc. Note that by doing so the condition of con¬ 
tinuity of slope of the displacement field (6) is not satisfied, 
but as we have seen, the solution is hardly affected by this. 

The second continuity condition (10) reduces now to 


V 

V, 


B 


ill 

*♦5 


(39) 


which means that | V I f- 1 V j . On the other hand formula (38) 

b L 

in the present solution implied that b f 0. Recalling that 
b = (V - V_) cosa, we have IV | = |V„', which contradicts the 
condition-of continuity. A consistent set of equations can 
be written, as discussed previously. However, the resulting 
solution, which apparently looses its appealing simplicity, 
is believed to introduce insignificant changes and thus will not 
be discussed further. 

The present solution with two moving hinges becomes increas¬ 
ingly inaccurate with decreasing value of the parameter Two 

factors are responsible for that. First, at an early stage of 
the crushing process, only one hinge circle is formed and the kine¬ 
matics of the process does not involve a free parameter to be opti¬ 
mized. Secondly, for small ot the portion of the toroidal shell 
becomes relatively shallow so that the assumption leading to 
l - 0 is no longer valid, and the shell should be treated as a 
three-dimensional body, according to the method developed by 
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Calladine [18]. Te shall study these two cases separately. 

Assume that the radius of the inner hinge is zero. From 
simple geometrical consideration Fig. 10/ it follows that 

b = R sina , R = 2 r 

Introducing (40) into (35) and (37) and eliminating R/h we 
obtain the following force-deflection characteristics 



4_ w 
3 h 


1 


which is identical to the one derived in [27]• The function 


(40) 


(41) 


(41) represents a straight line intersecting the curve described 
Ly (34) at g = |, Fig. 11. 

Return now to the general case of the shell in which the 


internal energy dissipation is governed by (13) and (16b). We 
assume that (40) holds, and that the shape of the shell middle 
surface represents a parabolic arc. The energy balance equation 
can be written in the form 


b 



■•here yl and yll are respectively the upper and lower contour of 
the shell cross-section, measured from the axis of instantaneous 
rotation, see Fig. 4. 



R 


(43a) 
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y 


ii 



(43b) 


The integrals in (42) can be readily evaluated and the 
crushing force becomes a function of a single parameter £ 


P 

2ttM 

o 




8 


15 R h 


Introduction (40) into (35) we obtain a unique relation 
between the radius of the hinge and central deflection 


,2 

w _ b 
h 2hR 

The minimization of the right hand side of (44) yields an 
optimum value of the parameter £ 


£ _1_ b^ 1 w 
h ~ 12 Rh 6 h 


(44) 


(45) 


(46) 


Substituting (45) and (46) back into (44) we finally obtain 


P 


2ttM 

o 


2 




( 47 ) 


The above relation is depicted in Fig. 11 by a solid line. 
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Ecuation (47) is valid as long as the line of instantaneous 
rotation first touches either the inner or cuter contour 
y^ cr v^ . This first occurs at ^ Thus, the 

range of applicability of both solutions appears to be the same. 
The formula (41) is easy to evaluate but it gives unsatisfactoy 
results for w/h = 0. This will become even more evident in the 
problem of a shell loaded ,’y a rigid boss. Cn the other hand, the 
solution (47) predicts correctly the yield point load but the 
mathematics involved are much more complicated. 

The 4% difference between the functions (34) and (47) at 
w/h = — is attributed to the fact that the former involves the 
simplifying assumption z, = 0, while the latter not. By con¬ 
tinuing the exact calculation beyond that point, the corresponding 
solution (doted line in Fig. 11) would gradually converge to the 
curve representing the formula (34). The thin full line in the 
same figure denotes the Morris and Calladine numerical solution. 
Unfortunately, no information was given in [17] on the magnitude 
of the integration step employed so it is difficult to assess 
its accuracy. A good agreement of both solutions is observed 
over the entire range of deflection covered by the Morris and 
Calladine solution. 

t. Crushing of a Spherical Shell Between Rigid Plates 

The procedure for solving this problem essentially follows 
that presented in the preceding section, except that now the 
rate of external work should be computed according to (17b) and 
the relation between w and a is different. As before, 
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the analysis follows several stages. The yield point and post¬ 
yield behavior is accurately predicted by the formula 


= 4 + 


! (h) 


w 3 

valid for ^ ^ . In the later stage of the crushing process, 

when the inner hinge leaves the shell axis, the plate displacement 
is related to the angle a by 


w R . 2 h 
_ = _ sin a _ _ sina 

The optimum value of width of the deforming radius is the same 
as before b = /3hR, Equation (38), and the resulting force- 
deflection relationship takes the form 

w _ _1_ / P . ,\ 2 _3_ / P . 

h “ 43 \2 ttM 0 ■ ) "32 \2ttM q 

The above solution is valid theoretically until the crush 

distance reaches the shell radius w = R. It is interesting to 

compare the formula (50) with an approximate solution due to 

Updike [14]: 



which was believed by the author to be valid for deflections 
not exceeding about one tenth of the shell radius. When both 
solutions are plotted, it is seen that a relative difference 
between them is small over the entire range of variation of w, 
Fig. 12. 


(49) 


(50) 


(51) 
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9. Rigid Boss Loading 

Leckie and Penny performed a series of tests on carefully 
manufactured spherical shells with R/’n = 60 and variable boss size 
[16]. They observed that with increasing downward deflection 
the crushing force for a given radius of the boss p reaches a 
maximum, then drops suddenly as the shell buckles, and starts 
to rise again. They also found that the initial plastic yield- 
point load depends markedly on a single geometrical parameter 

B = (52 

/Rh 

but with increasing deflection this dependence is much less 
pronounced. However, the tests were not run far enough to decide 
whether all P - w curves did converge to a single one. 

Assuming a rigid perfectly plastic material idealization 

Morris and Calladine [17] were able to predict accurately the 

initial response of the shell. Their numerical calculations have 

shown that while the slope of the load-deflection curves becomes 
w 

the same for ^ _> 2, the curves corresponding to different 
3 are considerably shifted with respect to one another. 

It is reasonable to expect that the effect of the boss 
size is a local one and should be almost entirely "forgotten" 
by the shell at later stages of the deformation process, so that 
all curves would eventually converge. Indeed, the deformation 
mechanism with expanding ring of plastic deformation is the same 
as before, except that the reference point for measuring deflec¬ 
tion changes with the boss size. This gives rise to a new term 
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in the force-displacement relationship 

2 

+ i: 


X = J_ /_L_ + A - 3 . 
h 16 \2 ttM + 2 2 


(53) 


The above solution reduces to the one describing the shell 
under central-point load (formula (34)) by setting £ = 0. The 
new term is relatively small, and provides a shift of P - w 
curves in the w/h direction Fig. 13. The horizontal displacement of 
P-v; curves in the numerical solution reported in [17] for the same 
values of 8 are much larger, which is probably attributed to 
the trial and error optimization procedure employed. 

We were unable to derive a simple closed-form solution for 
the initial phase of loading based on the exact formula (13) 
showing the nature of the convergence. However, the interesting 
pattern of behavidr of the shell associated with plastic snap- 
through can be described, at least qualitatively using the 
approximate formula (14). Assuming that a stationary hinge circle 
is formed at T = p and the moving one at T = p + b, the load- 
deflection relationship was found to be 


P 

2itM 



-B 


Ft* 


2 B 


FT 


+ 6-8 


1 


(54) 


The above solution reduces to (41) for 8=0 but it leads to an 
unrealistic infinite force magnitude with ^ 0. This is a 

consequence of the assumed symmetric geometrical shape of the 
shell cross-section (36) with zero initial width of the plastic 
zone to initiate the motion. In reality, the shape is unsymmetric 
with finite width. It was possible to derive an approximate 
expression for the initial yield-point load (w/h = 0) in the form: 
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2ttM 


= 86, b = 


/Rh 


Rh 

2p 


(55) 


The plot of the function (54) with (55) for several values of 
the parameter 6 is shown in Fig. 14. 

10. Spherical Cap Under External Pressure Loading 

Consider a simply supported spherical shell with base radius 
a, and rise H, Fig. 15. The shallowness of the shell is described 
by the parameter 


a 



(56) 


Following Ref. [21] define the total load P and the reference 
load P Q by 

P = -rpa 2 P = 6ttM (57) 

c o o 

where p is a uniformly distributed downward pressure. It is 
relatively easy to derive the load-deflection relationship for 
any value of a, using the present method. We shall however 
restrict the analysis to the limiting cases of small and large 
a. Consider first the case a >> 1 for which equation (14) applies. 
Suppose the outer boundary of the plastically deforming zone BC 
is b while the portion AB remains rigid. The assumed velocity 
field consist, as previously, of the rotation around the point B 


• • 

w(x) = tux (x 1 b £ a ) (58) 

This rotation imposes a continuous change in the initially 

circular shape of the shell middle surface. If the central 

deflection is smaller than the rise of the shell above the 
b 2 

circle b, w < , the current shape can be approximated by 
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the parabolic arc. 

y(x) = (l)~ (Ir + w ) 

For the time being we assume the plastic hinge circle to be 
located somewhere between the shell center and the support 
0 < b <_ a. Equating the rate of external work to the energy- 
dissipated in the region of continuous deformation and plastic 
hinge circle one gets 



2 7T 


/ 


pw(x)xdx = 2irq 


7 0 ha) / y (x 


) dx + 2tto q — mb 


or after integration 


16 0 ) 


P_ 

P 

o 



i Si + i 

3 h + 1 


(61) 


Id 

where b = JhR' We can see for t * ie force to be a minimum 

one has to consider the largest possible value of b, i.e. b = a. 

With the assumed simply supported boundary conditions there will 
not be any contribution of the concentrated hinge to the total 
dissipation energy, and the last term in Eq. (61) should vanish, 
giving the final expression for the force deflection characteristics 

— = ja - for w < a (6 2) 

o 

In particular, the initial yield-point load is obtained from 
(62) by putting w = 0, P/P 0 = 8/3a. A plot of the above formula 
is undistinguishable from the curve representing the much more 
complicated solution derived by Kondo and Pian [21], Fig. 16 

At the other extreme we have very shallow shells in which 
the line of instantaneous rotation lies within the contour of 
the shell cross-section, Fig. 4. This implies a < 1. 
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Assume that the shape of the shell middle surface is 
changing with deflections according to Eq. (59). The energy 
in the continuous deformation field is given by 


E , = 2ira_ 

cont. o 


f h \ r yI / 

~y ZI ) 

/ J J 

ydy j 


where the limits of the integrals appearing in (63) are 

I h r 

y = Y + 2 “ ? 


II ^ h ^ r 

y = -y + , + C 


The expression for the energy dissipated at the outer hinge 
would depend on the position of the hinge circle. If b < a, 
then the boundary conditions do not have any effect on the 
solution and Eq. (16b) applies. If the plastic zone extends 

up to the support, b = a, then for a simply supported shell 
restrained from axial motion, the term responsible for bending 
vanishes 


E, . = 2TTG'_tok 

hinge c 


r h 2 + 


n = 


J 


1 if b < a 
(0 if b = a 

The left hand side of the energy balance equation is the same 
as before (60), while the right hand side is a sum of (63) and 
(65). Evaluating the integral, one obtains 


P b 

P a 
o 


(1 + H) + ||b 2 -^bq + |-q 2 + 2C 2 - : j5(2b-q) 


where 



* q = 


w 

h 


(63) 

(64a) 

(64b) 


(65) 


( 66 ) 
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After minimization with respect to £, the formula ( 66 ) becomes 


P b , 56 r2 68 t— 26 2 

p a ~ (1 + n) + 45 b 45 bq 45 q 

o 

The load-deflection curve, represented by (67) is dropping 
with increasing b and attains the lowest position when 


(67) 


-2 45 A 26 , 

b = (1 + n) 56 + 5 ^ g 2 


( 68 ) 


For example at q = 0, the optimum value of radius of the outer 
hinge is b = 1.3, falling beyond the range of applicability of 
the present solution, which is b < a < 1. Consequently, the least 
value of the crushing force is obtained by taking the largest 
possible radius, i.e. b = a. We have thus shown that for shallow 
shells the zone of continuous plastic deformations is extending 
right up to the support. Substituting n = 0 and a = b into (59), 
leads to the following force-deflection relation. 


_P 

P 


56 2 
45 


— = 1 + — OL - Tr- Ctq + — q 


68 

45 


26 _2 
45 


(69) 


In particular, the initial yield-point load is given by 

P , . 56 2 . 

P“ =i 45 a a<l 

o 


(70) 


In deriving Eq. (69), shell deformations both in radial and 
circumferential directions were taken into account. A relative 
contribution of the two components of the dissipated energy is 
controlled by the position of the axis of instantaneous rotation I~I. 
The optimum position is one which minimizes the magnitude of 
the pressure required to maintain the plastic flow. It is 
interesting to note that the neglect of meridional strain, an 
assumption made in Ref. [21], would lead tc an increase rather 
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than decrease in the level of the yield-point load. Indeed, in 
the present approach all work done by the meridional strain is 
concentrated in the generalized tensile hinges and the corres¬ 
ponding dissipation energy is given by (65). This term can 
easily be made equal to zero by taking 5=0. Substituting 
a = b and 5=0 into ( 66 ) we obtain a new solution 


P 

P 


. x 32 r2 12 ,4 2 

= 1 + — b _ -F- aq+^-q 


(71) 


o - 15 “5 

Load-deflection curves for few chosen values of the parameter a 
are shown in Fig.17 together with the solution due to Kondo and Pian. 
The initial yield-point load is obtained from (71) by setting 


q = 0 : 


p i 32 2 ^ 

= 1 + rj-jr a a £ 1 <72) 

o 

The above formula is identical to the one derived by Kondo 
and Pian and the corresponding P - a curve clearly lies above 
that given by (70), Fig. 16. Thus, the neglect of meridional 
strain rate was offset by even larger increase of the dissipa¬ 
tion due to the circumferential strain rate. With increasing 
ct, the contribution of dissipated energy in meridional direction 
diminishes, illustrating a nature of the approximation introduced 
in replacing (13) by (14). Note that the solution (62) was derived 
on the assumption that 5 = and this explains the excellent 
correlation of the present solution with the more elaborate 
solution due to Kondo and Pian, which employs a full set of equations 
describing moderately large deflection of shells. 









-34- 


It is not in principle difficult to extend the present 
solution to shells with intermediate values of the parameter a. 
The P - a curve would first follow Eq. (70) and then converge 
assymptotically to Eq. (62). 


11. Effect of End Condition 

According to the present analysis, the boundary conditions 
are not felt by the shell until the outer hinge reaches the 
support located at p*. It is seen from (35) and (38) that this 
occurs when deflections become 



(73) 


For a shell under•central-point load, the width of the plastic 

zone predicted by the approximate solution is b = /3Rh. From now 

on the outer hinge remains fixed, and any further increase of w/h 

is due to the diminishing of the width of the toroidal section b. 

Introducing the dimensionless parameter a = —^ ■ , the geometrical 

/3Rh 

relation becomes 


w_p* 2 32 w* 3 

h ~ TT ~ 2 a = TT + 7 


(1 


2 


a 


) 


For fully clamped boundary conditions, the solution (37) with 
two concentrated plastic hinges still applies, so that after 
introducing (73) into (37) we have 


P 

2ttM 0 



/3Rh 


+ b 


i 


(74) 


(75) 


Equations (74) and (75) provide a parametric representation of the 
load-deflection relationship at the terminal stage of the shell motion 
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Kith a -*■ C the force blows up to infinity .Fig- 18. This occurs 
when + j. In reality the force level in shells 

subjected to central point load is equal to the strength of the 
inverted spherical shell in tension. An upper bound for this 
force is 


P tens _ 4R 
2 M q t h 


(76) 


This is much higher than (75) with a - 1, but finite. For shells 
compressed between rigid plates the formula (75) applies 
multiplied by the coefficient 2. The force increases with a 
approaching zero and indeed goes to infinity when the upper plate 
comes in contact with the supporting plate. A typical loading 
situation for the two cases described above is shown in Fig. 18. 

12. Conclusions 

This paper presents a general methodology for predicting the 
crush resistance of arbitrary rigid plastic shells of revolution 
deforming in an axisymmetric mode. The implications of the condi¬ 
tions of kinematic continuity, as well as the boundary conditions, 
are discussed. Several practical problems are solved, including 
a conical shell and a spherical shell crushed between rigid plates 
and under boss loading, and a spherical cap under external uniform 
pressure. 

For the initial stage of deformation the structure has been 
treated as a three dimensional body, in the spirit of Calladine's 
approach [18], while in subsequent stages a simplification to this 
approach has been suggested, which results in simple closed form 
solutions useful for engineering applications. A good correlation 
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was found with some existing solutions of similar problems obtained 
through much lengthier calculations. 

In all cases treated, the zone of plastic deformation was 
found to be confined to a relatively narrow region. Also, the 
level of the crushing force was shown to be quite insensitive to 
the particular shape of that 2 one. 

The present method makes it possible to study the crush 
behavior of rotationally symmetric shells strengthened by circum¬ 
ferential ring stiffeners. The stiffner acts initially as a 
clamping ring until the force reaches a critical level. Then the 
ring becomes a rigid boss. Thus the determination of the initial- 
point load for a shell loaded by a rigid boss provides a basis 
for studying the optimum design of stiffened plastic shells 
* against crash. In the optimum design the rigidity of the stiffeners 
must be compatible with the local crushing strength of the shell. 
This problem will be studied in detail in a future publication. 
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